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WALSH-FUNCTT ON  REPRESENTATION  OF  SHIFT  REGISTERS 
WITH  STOCHASTIC  INPUTS 

A. Udaya  Shankar  and  David  K.  Cheng 
Electrical  and  Computer  Engineering  Department,  Syracuse  University, 

Syracuse,  New  York 

ABSTRACT 

A WaJ sh-Hadamard  description  for  the  operation  of  linear-feedback 
shift  registers  is  presented.  It  is  used  to  study  the  transient  and 
steady-state  behavior  at  the  outputs  of  shift  registers  with  respect  to 
a stochastic  input  and  an  initial  register  state. 

Introduction 

Despite  the  wide-spread  use  of  linear  feedback  shift  registers 
(LFSRs) , a method  for  determining  the  transient  and  steady-state  output 
statistics  of  an  LFSR  in  a noisy  environment  under  an  initial  register 
state  does  not  appear  to  be  available.  In  this  article  we  make  use  of  a 
Walsh-Hadamard  representation  for  the  outputs  and  the  external  input  and 
characterize  the  LFSR  operation  by  a set  of  time-recursive  equations  over 
the  reals.  The  equations  can  be  arranged  in  a matrix  form  where  the  matrix 
involved  is  sparse  and  can  be  written  by  an  inspection  of  the  LFSR.  Given 
a description  of  the  LFSR's  initial  state  and  external  inputs  that  are 
white  and  stationary,  the  Walsh-IIadamard  representation  allows  an  easy 
computation  of  the  transient  and  steady-state  output  statistics  and  dis- 
plays the  convergence  of  the  transients  in  an  elegant  manner. 

PreJ  imi.nar ies 

Figure  1 is  a schematic  of  an  n-stage  LFSR  , where  the  feedback 
coefficients  a_,  a.,...,  a . are  0 or  1 with  a 1 indicating  a tap.  To 


1 w 


ensure  a genuine  n-stage  LFSR  we  insist  that  a , = 1.  The  feedback  at 
i n— 1 

n-1 

time  t is  ^ a y (t-1),  where  denotes  modulo-2  addition.  The  transi- 
k=0  k k 

tions  of  the  LFSR  are  then  described  by 


n-1 

y (t)  = x(t-l)  © § ay  (t-1) 

° k=0  k k 


yR(t)  = yk_1(t-l),  for  0 < k < n 


(1) 


(2) 


For  any  variable,  say  x,  ranging  over  {0, 1 , . . . , 2n-l } , we  denote  its 

binary  representation  <x  x,,x>  by  x.  On  the  domain  of  the  binary 

n-1  1 0 

2 

n-tuples,  the  ith  Walsh-Hadamard  function  is  defined  as 


n-1  n-1 

A HH  / Hh 

Wal[i,xJ  = (-I)*"  0 = (-1)*  0 


for  0 < i,  x < 2n 


(3) 


n-1  n-1  , 

where  i = l ijj,2^  and  x = £ x„2 

£=0  1=0  * 


Because  the  set  of  all  Walsh 


functions  on  the  n-tuples  form  a complete  orthogonal  basis,  any  problem 
on  the  binary  n-tuples  can  be  solved  equivalently  in  the  corresponding 
Walsh  domain. 


Walsh-Hadamard  Description  of  LFSR 
Using  the  quantities  in  eqns.  1 and  2 as  the  exponents  of  (-1), 


we  obtain,  respectively 


In  view  of  eqn.  3,  eqns.  4 and  5 become 


Wal[2k,  y(t) ] = / 


Wal[l,  x(t-l)]  Wal[a,  y(t-l),  for  k = 0 


Wal [ 2k_1 , y(t-l)],  for  0 < k < n 


(6) 

(7) 


where  y(t)  = <y  . (t) , . . . ,y, (t) , y„(t)>  and  a = <a  , ...,a  ,a  >.  We 
n— i l u n-i  1 U 

note  that  cqns.  6 and  7 are  insufficient  to  sustain  recursion  because 
Wal[a,y(t)]  is  not  generated. 

Now,  for  every  subset  of  eqns.  6 and  7 a product  can  be  formed  to 
yield  a unique  equation  relating  y(t)  and  y(t-l).  Let  i for  1 i < 2n, 


indicate  by  its  l's  the  equations  from  eqns.  6 and  7 which  are  multiplied. 
Equating  the  resulting  left-hand  and  right-hand  sides,  we  have 


II  {Wal  [2  y (t)  ] } = {Wal  [ 1 ,x(  t-1)  ] Wal[a,y(t-1)  ] } ° 

«,= 0 


x n {Waif 2^  1,  y (t-1) ] } £ (8) 

£=1 


il 


For  convenience,  we  denote  Wal[i,  y(t)]  by  Y^(t)  and  Wal[l,  x(t)]  by 
X(t).  From  the  property  that  Wal[i,  y(t)]  Wal[j,  y(t)]  = Wal[i  © j»y(t)], 
where  © is  extended  to  denote  dyadic  addition  (component-wise  modulo-2 
addition  on  i and  j),eqns.  8 yield 


4 


[Y^t),  Y^(t) YN(t)]T  by  Y(r),  where  N - 2U-1 , we  can 

9 and  10  In  a matrix  form 


ins . 


Y(t)  = B [X(t-l) ] Y(t-l)  (11) 

where  the  N * N ’-matrix  B[X(t-l)]  = {B_  [X(t-1 ) ] : 1 < i,  j £ N)  has  the 

following  elements: 


1,  for  even  i and  j = i/2 


B . [X]  = ( X,  for  odd  i and  j = (-“^)  © a 


(12) 


,n-l 


2°  ^ < (— —•)  © a _<  N since  ^ 


0,  otherwise 

For  1 _£  i N:  if  i is  even,  then  1 i/2  < 2“  and  if  i is  odd, 

= 1. 

It  is  clear  from  eqn.  12  that  B[X],  for  X / 0,  is  a sparse  matrix 
with  exactly  one  nonzero  entry  in  each  column  and  in  each  row  and  that 
the  entry  is  1 for  the  first  2n  ^-1  columns  and  X for  the  2n  ^ succeeding 
columns.  Thus,  B[X]  is  essentially  a permutation  matrix  involving  the 
scaling  by  X of  some  elements.  For  example,  a 3-stage  LFSR  with  feedback 
coefficient  vector  a = <1  0 1>  will  have  a 7 x 7 B-matrix.  Since  1 © a = 4, 


a = 7 

and  3 

© 

a = 

6, 

we 

have  from  eqn. 

11,  using 

eqn.  12, 

Yj(t) 

- 

0 

0 

0 

0 

X(t- 

-1)  0 

“S 

0 

Y1(t-1) 

Y2(t) 

1 

0 

0 

0 

0 

0 

0 

Y2(t-1) 

Y3(t) 

0 

0 

0 

X(t- 

-1)  0 

0 

0 

Y.(t-l) 

Y4(t) 

0 

1 

0 

0 

0 

0 

0 

Y. (t-1) 

4 

Y (t) 

0 

0 

0 

0 

0 

0 

X(t-l) 

Y5(t-1) 

0 

0 

1 

0 

0 

0 

0 

Y6(t-1) 

*7(t) 

0 

_ 

0 

0 

0 

0 

x(t-l) 

0 

Y?(t-1) 

(13) 


Returning  to  eqn.  11,  we  can  express  Y(t)  in  terms  of  its  original 


state,  Y(0) . 


I 


Y(t)  = B[X(t-l)]  B [X(t-2)  ] . . . B [ X(0) ] Y(0)  (14) 

The  product  B[X(t-l)]  B[X(t-2)]...  B[X(0)]  is  clearly  also  a permu- 
tation matrix  having  1 or  -1  as  its  nonzero  entries.  When  x(t)  = 0, 

implying  X(t)  = 1,  the  LFSR  operates  as  an  oscillator  with  a sequence 

e 

period  which  we  will  denote  by  p.  In  this  case,  y(p)  = y(0)  and  hence 
B^[l]  = I where  I is  the  identity  matrix  of  an  appropriate  dimension. 

Output  .Statistics 

We  assume  a white  stationary  input  so  that  x(t^)  is  independent 
of  but  identically  distributed  as  for  t^  ^ t2*  This  allows  the 

description  of  the  initial  state  of  the  LFSR  to  be  independent  of  the 
input  x(t)  for  t ^ 0.  With  E [ * ] denoting  expectation,  let  ijj^(t)  stand 
for  E[Y^(t)],  vjj ( t ) for  E[Y(t)],  and  y for  E[X(t)].  Then,  taking  the 
expectation  of  eqn.  14  results  in  the  following  transient  solution. 

$»(t)  = Bt[X]ip<0)  (15) 

The  steady-state  statistics,  if  they  exist  (i.e.,  if  Bt[x]  converges), 
can  be  found  by  equating  ijj(t)  to  ijj(t-l)  in  the  expectation  of  eqn.  11. 

We  have 

^(t)  (I  - B [ xl ) = 0 (16) 

It  is  seen  from  eqn.  16  that  if  (I  - B[xl)  ' exists,  the  steady-state 
ijj(t)  is  uniquely  0.  This  means  that  al]  the  outputs  are  independent 
and  equiprobably  0 or  1.  Obviously,  for  the  oscillator  case  where 
X = 1 or  -1,  a steady-state  statistic  does  not  exist  except  for  the 
trivial  case  of  iKO)  = 0.  Thus,  (I  - B[l])  and  (I  - B[-l])  do  not 


exist. 


The  eigenvalues  of  Ii [ x 1 control  the  time  u-igenet  . tin  LFSK 


A 


output  statistics.  Indeed,  denoting  the  N eigenvalues  of  B[xl  by 
* * " *^N*  some  wbich  may  be  equal,  we  may  write 

B^xl  = Aj  C(l)  + A*  C(2)  +...4  A*  C(N)  (17) 

where  the  matrices  C(l),  C(2),...,  C(N)  are  determined  from  the  eigen- 
vectors. As  t -*•  an  eigenvalue  of  a magnitude  greater  than  1 would 
cause  Bt[x]  and  hence  tp ( t ) to  explode,  which  is  not  allowed  physically. 

Thus  |X  | <_  1 for  1 £ i < H.  If  the  magnitude  of  an  eigenvalue  is  1, 
eqns.  15  and  17  indicate  that,  at  steady  state,  $(t)  would  be  a non- 
decaying periodic  function.  In  view  of  the  constantly  growing  uncer- 
tainty introduced  by  x(t),  such  a ijj(t)  is  possible  if  and  only  if 
X E {1,  -1).  This  implies  determinism  right  from  the  outset.  Thus,  if 

|x|  < 1,  we  expect  the  contents  of  LFSR  to  ultimately  reach  maximum 

entropy;  i.e.,  ijj(t)  H 0. 

The  above  observations  may  be  summarized  as  follows:  (i)  | A_J  £ 1 

for  all  1 i < N;  (ii)  | x I = 1 if  and  only  if  |A_J  e {1,0}  for  any 

1 < i < N with  at  least  one  eigenvalue  having  a unit  magnitude;  and 
(iii)  |x|  < 1 if  and  only  if  | A^|  < 1 for  all  1 <_  i < N.  From  eqn.  17, 
if  A is  the  eigenvalue  of  the  greatest  magnitude  and  |x|  < 1,  then  the 
entries  in  Bt[xl  will  decay  to  zero  at  least  as  fast  as  Um|t- 

The  joint  probability  distribution  of  y(t)  can  be  obtained  from 
{}}( t)  via  the  Walsh-transf orm  relationship.  By  definition, 

N 

iMt)  = l Pr[y(t)  = l]  WalfiJ]  (18) 

1 £=0 

The  inverse  transformation  of  eqn.  18  is 


11  08  015  9 


w 


Pr[y(c)  = M]  = 2_n[l  4 l ip  (t)  Wal [ i,  |"I]  ] (19) 

i=.l 


We  may  conclude  that,  for  | x I < 1»  C t)  converges  to  0 and  Pr[y(t)  = jj] 

converges  to  2 n,  also  at  least  as  fast  as  |X  I . 

' 

For  the  example  considered  in  eqn.  13  the  output  statistic  Y(t)  = 

Pr [y(t)=<001.>]  is  plotted  versus  t in  Fig.  2 for  an  initial  state  <001> 
two  values  of  y*  For  this  maximum-length  3-stage  LFSR  (p=7) , = 

X4/?  exp(2Tri  /-1/7),  i = 1,2,..., 7.  Thus,  | A ^ j = X ^ ■ We  observe  from 
Fig.  2 that  (i)  Y(t)  = 0 for  t = 1 and  t = 2 irrespective  of  x(t)  for  the 

_3 

given  initial  state  <0  0 1>,  as  expected,  (ii)  y(t)  converges  to  2 = 

0. 125,  as  predicted  by  eqn.  19;  and  (iii)  |y(t+7)  - 2 ^ | / |y(t)  - 2 = x* 

Convergence  is  faster  for  a smaller  y. 

Conclusion 

The  behaviour  of  linear-feedback  shift  registers  is  neatly  captured 
by  a Walsh-Hadamard  representation  in  the  form  of  a matrix  recursive  equa- 
tion which  can  be  written  by  inspection  and  which  yields  the  output  sta- 
tistics for  a stochastic  input  quite  easily.  Because  the  representation 
is  over  the  reals  where  our  intuition  is  stronger,  it  facilitates  our 
study  of  the  transient  and  convergence  properties  in  terms  of  the  eigen- 
values of  a permutation  matrix. 
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Fig.  2 - Output  probability  of  a 3-stagc  maximum-length  LFSR 

a = <1  0 1>,  y (0)  = (i  = <0  0 1>, 


versus  time: 
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